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ABSTRACT 

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium 
(ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS 
spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [01] 63 
yum, [OIII] 88 fim, [Nil] 122 yum, [CII] 158 //m and [Nil] 205 fim. We observe that in the [OI] 63 yum, [OIII] 88 yum, and [Nil] 122 fim line maps, 
the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation 
velocity of ~ 220 km s"' (inclination uncorrected) measured in all lines. The [OI] 63 pm/[CII] 158 yum ratio varies smoothly throughout the central 
region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks 
in the [01] 63 yum and [OIII] 88 fim maps. Variations of the [Nil] 122 //m/[NII] 205 fim ratio correspond to a range in the ionized gas density 
between 150 and 400 cm"^^. 
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[ 1. Introduction 

NGC 1097 is a Seyfert 1 galaxy with a bright starbu rst ring with 

a diameter of 2 k pc and a strong large-sc ale bar (J Gerin et aO 

i |l988HKohnoetal. !2003: H sieh et al.ll2008h with a length of 15 

' kpc. Optical and near-infrared images reveal dust lanes that run 

along the primary large-scale (15 kpc) bar and curve into the 

ring, which is form ed by two very tight spiral arms, and a second 

_ bar inside the ring dOuillen et al. 1995; Prieto et al. 2005) . This 

, bar may be responsible for driving gas into the n ucleus, pos- 

■ sibly fueling the central su per-massi ve black hole dPrieto et al.l 

; [2005; Fathi et al. 2006; Dav ies et al.H 2009), and may also have 

triggered the formation of a compact star cluster seen near the 

nucleus. NGC 1097 provides an excellent opportunity to study 

the physical conditions of the interstellar medium (ISM) in a 

. nearby galaxy with both a starburst and an active nucleus. 

The nucleus and the star- forming ring are prominent in CO 



the ring, along with faint line emission from the inner region, 
exhibiting a LINER-type spectrum. 



ling ring a 
and HCN line emission (|Kohno et al.ll2003h. Near-infrar ed spec 
troscopy (iReunanen et alJl2002t iKotilainen et ani20 00^ reveals 



emission from both ro- vibrational H2 and H-recombination lines 
at the nucleus, the star-forming ring, and the region in between . 
Optical long-sUt spectroscopy (iStorchi-Bergmann et al.l [l996l) 
shows strong ionized gas emission from both the nucleus and 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
** Hubble Fellow 



With Herschel/PACS ("Poglits ch eraLll2010t) we are now able 
to target the most important cooling lines of the warm ISM on 
physical scales much smaller than ever before possible in exter- 
nal galaxies. The KINGFISH project (Key Insights on Nearby 
Galaxies: a Far- Infrared Survey with Herschel - PI: R. C. 
Kennicutt) is an open-time Herschel key program which aims 
to measure the heating and cooling of the gaseous and dust com- 
ponents of the ISM in a sample of 61 nearby galaxies with the 
PACS and SPIRE instruments. The far infrared spectral range 
covered by PACS includes several of the most important cool- 
ing lines of the atomic gas, notably [CII] 158 yum, [OI] 63 //m, 
[OIII] 88 yum, [Nil] 122 jum, and [Nil] 205 /urn. NGC 1097 is 
one of the KINGFISH targets selected for the Herschel Science 
Demonstration Program (SDP) (for PACS imaging of NGC 1097 
see Sandstrom et al. (2010) in this volume and for SPIRE obser- 
vations see Engelb racht et alJ (120101) in this volume). 

In this letter, we present far-infrared spectral line maps 
of the circumnuclear starburst ring and the large-scale bar of 
NGC 1097, obtained with the PACS spectrometer on board the 
ESA Herschel Space Telescope (Pilbratt et al. 2010). The maps 
presented in this letter are the first PACS spectral maps from 
the KINGFISH program. Througho ut this paper we as sume a 
distance to NGC 1097 of 19.1 Mpc dWilhck et alJ[T997l) . which 
gives a projected scale of 1" - 92 pc. 
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Fig. 1. Composite of [CII] 158 /iin emission maps in the 
center (yellow contours) and extranuclear positions (red con- 
tours), overlaid on a Spitzer/MlPS 24/7m image of NGC 1097 
(in grayscale). The extranuclear maps are a combination of 
wavelength-switching data (North) and chop-nod data (South). 
The contour levels in units of 10"^ W/m^/sr are: 0.35, 0.70, 1.4, 
2.8, 4.2, and 5.95 for the nucleus, 0.092, 0.184, 0.276, 0.368, 
0.460, and 0.570 for the southern extranuclear region, and 0.05 1 , 
0.102, 0.153, 0.204, and 0.238 for the northern exti-anuclear re- 
gion. 



2. Observations and data reduction 

NGC 1097 was observed with the 5x5 pixel integral field 
unit (IFU) of the PACS Spectrometer in both chop-nod (CN) 
and wavelength switching (WS) modes as part of a series of 
Science Demonstration Phas e tests. The field of vie w of the IFU 
is 47" X 47" with 9'.'4 pixels (iPoghtschetal. 11201 Oh . The diffi-ac- 
tion Hmited FWHM beam size of PACS is 5'.'2 at 60 - 85//m 
and 12" at 130 - llOfim. The spectral resolution of PACS is 
180 kms-i for [OI] 63 yum, 120 kms' for [OIII] 88 yum, 290 
kms-' for [Nil] 122 /urn, 240 kms"' for [CII] 158 jjm, and 150 
kms"' for [Nil] 205 fim. CN observations were carried out be- 
tween January 4 and January 22, 2010, with an observing time of 
19.4 hrs, with 17.7 hrs on-source. WS observations were carried 
out between January 4 and January 22, 2010, and the observing 
time was 6.6 hrs, of which 3.6 hrs are on-source. Five different 
regions were observed: the nucleus, two extranuclear positions 
at the outer ends of the bar, and one strip at each side of the nu- 
cleus, following the minor axis of the galaxy. CN observations 
were obtained in a 2 x 2 raster in the nucleus for all five target 
lines (with a step size of 9'.' 4), in a single short strip along the 
minor axis ([OI] 63 yum, [Nil] 122 //m, and [CII] 158 //m), and 
in a single extranuclear position to the south ([OI] 63 /zm, [OIII] 
88 yum, [Nil] 122 /zm, and [CII] 158 yum). Each map covers an 
area of 4.5 x4. 5 kpc^. In the basic CN mode, the telescope chops 
between two regions of sky 6' apart, with the target placed alter- 
nately in one or the other beam. WS observations were obtained 
for the same lines in a 2 x 2 raster in the nucleus, two short 
strips along the minor axis on each side of the nucleus, and in 
two extranuclear positions to the north and south of the nucleus. 
Line fluxes obtained from WS data in the same regions were in 
all cases comparable within the flux uncertainties to CN results, 
and the KINGFISH survey will utilize WS for the remainder of 
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Fig. 2. Spectral maps of the lines observed with PACS at the nu- 
clear position, and an image of the ring at 70 yum with PACS. 
The maps represent the integrated flux over -500 <Vsys < 500 
kms"', with v^s - 1271 km s"'. North is up and East is to the 
left. In all cases the size of the image is 51" x 51". The cross 
marks the location of the radio position of the nucleus. 



the sample. Here we focus on the CN observations, which have 
the best signal to noise ratios. 

The spectra were reduced using HIPE {Herschel Interactive 
Processing Environment), Special Development version 3.0 CIB 
1134. The pipelines remove detector artifacts, obvious cosmic 
ray signatures and apply ground-based flat-field corrections. 
Dark-current subtraction was done for the WS case. Additional 
band-dependent scale factors were applied to the post-processed 
extracted fluxes to adjust the ground-based absolute calibration 
to in-flight observations of calibrator^ The uncertainty in the 
line fluxes is dominated by the absolute flux calibration. The flux 
calibration uncertainties are on the order of 30%, and pixel-to- 
pixel relative calibration uncertainties are on the order of 10%. 

3. Results 

In Fig. 1 we present an overlay of the [CII] 158 yum maps of 
the nuclear and extranuclear regions of NGC 1097 on a 24 yum 
image taken with the Mu ltiband Imaging Photometer for Spitzer 
MIPS (iRieke et al.ll2004 ). The 24yL(m emission is well traced by 
the [CII] 158 pva emission in all the regions where [CII] 158 
jim emission is observed. As there are no CN observations of 



' Spectral extractions followed the prescription described in the 
Herschel PACS Data Reduction Guide and post-launch corrections to 
the "nominal response" followed the guidelines of Version 1.0 of the 
Spectroscopy and Calibration PACS release note. 
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Fig. 3. Integrated spectra of the nucleus and ring. The system 
velocity Vj,.j = 1271 kms"' is set as the zero-point, marked by 
the dashed-dotted line. The dashed line is the spectrum of the 
[CII] emission in wavelength-switching mode. 



the northern extranuclear region, we show a WS map, which 
clearly follows the detailed structure evident in the 24fim image. 
The signal-to-noise ratio in the minor axis strips is very low, and 
these data are not shown in Fig. 1 . 

In Fig. 2 we present continuum-subtracted maps of the cen- 
ter of NGC 1097, for the [CII] 158 fim, [OI] 63 jum, [OIII] 
88 jj.ni, [Nil] 122 //m, and [Nil] 205 jum emission lines. For 
comparison, we also show a PACS 70/vm continuum image (see 
ISandstromet"aD (12010) in this issue). The [OI] 63 i^m, [OIII] 88 
yum, and [Nil] 122 yum maps are the most similar There is a par- 
tial clumpy ring-like structure, with peaks NW, N and E of the 
nucleus in the [OI] 63 fim and [OIII] 88 fim maps, whereas the N 
clump is absent in the [Nil] 122 /im map. The peaks observed in 
the [CII] 158 fim coincide with the clumps observed in the [Nil] 
122 jjm map, but the overall distribution is much smoother. Only 
the N and NW peaks have a counterpart in the PACS lOj-im im- 
age. The [Nil] 205 //m looks strikingly different than the other 
maps. It has a peak NE of the nucleus but lacks the NW and E 
hotspots present in the other maps. The lack of a resolved ring 
in the [CII] 158 yum and [Nil] 205 /im may be partially due to 
the large beam at these wavelengths. A careful measurement of 
the spectrometer PSF on 10" scales will be required in order to 
further separate real physical components in a system as small 
as the NGC 1097 ring. 

To extract integrated spectra from the line maps, we selected 
a circular region of 23'.'8 in radius centered on the nucleus which 
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Fig, 4, Velocity field of the [OI] line in the center of NGC 1097. 
The contours represent approaching (blue) and receding (red) 
emission relative to the systemic velocity, \s\-s = 1271 kms"'. 

Table 1. Integrated line fluxes' (10"'^ Wm"^) 

Position [OI] [OIII] [Nil] [CII] [Nil] 

63 fim 88/im 122yum 158/jm 205/im 

0.15 



Nucleus 


3.5 


0.7 


2.0 


9.0 


Enuc S 


0.28 


0.13 


0.07 


1.2 


Enuc N 


- 


- 


- 


0.9 



'"Nucleus" comprises a circular region of 23'.'8 in radius (~ 2 kpc) 

centered on the Seyfert nucleus. The fluxes for the other positions were 

measured over the whole maps (50" x 50" ~ 4 x 4 kpc). All measured 

fluxes have an uncertainty of 30%. 



includes the ring, but avoids the noisy edges of the map. The 
continuum subtracted velocity profiles for each line are shown 
in Fig. 3. In all cases, zero velocity was defined as v,,., = 1271 
kms '. An overplot of the [CII] 158 //m velocity profile made 
from the WS data is also shown for comparison, which has a 
peak flux of ~ 90% of the CN peak, within the flux uncertainty. 
The smooth appearance of the [CII] 158 /im and the [Nil] 122 
/im line profiles is due to the worse spectral resolution at these 
wavelengths. The small difference is attributed to transients in 
the CN data which are not present in the WS data. The most 
notable characteristic of these profiles is that they are all double- 
peaked, which is expected for a rotating ring. In Fig. 4 we show 
the velocity map of [OI] 63 //m. The velocity spans a range of 
±220 kms ', being redshifted to the SE, and blueshifted to the 
NW. The [OI] velocity field is consistent with circular rotation 
also seen i n ionized (Ha, Fat hi et al. (2006)) and molecular gas 
(CO (2-1) iHsiehet all (120081) ). In Table 1 we present the inte- 
grated line fluxes for the nucleus and extranuclear positions. 

4. Diagnosing the ionized and neutral gas 
properties 

With the line maps we can study spatial variations in the prop- 
erties of the warm atomic gas. Fine structure line ratios trace 
the intensity of the incident radiation field Go on the neutral gas 
and the electron gas density n^ , respectively (e.g. jKaufman et aP 
1999). In Fig. 5 we show maps of the [OI] 63 /im/[CII] 158 /im 
and the [Nil] 122 //m/[NII] 205 //m line ratios, overlaid with 
the contours of the PACS 70/im image. Before making the ratio 
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maps, each image was smoothed to the beam size at the longer 
wavelength using a Gaussian profile, while conserving flux. The 
[Ol] 63 yL(m/[Cll] 158 fim ratio varies smoothly with values rang- 
ing between 0.25 - 0.45 throughout the central region, and is en- 
hanced on the northeast part of the ring. This enhancement par- 
tially coincides with the peaks in the [Ol] 63 yu m and [OIII1 88 
^m li ne maps and the peak of the Ho- emission (Humm el et al.l 
119871) . which indicates that the [Ol] 63 //m/[Cll] 158 //m ratio 
traces the most massive star forming knots in the ring. The values 
of [Ol] 63 /im/[Cll] 1 58 um are a factor of 2 lower than for the 
starburst galaxy M82 (IColbert et al.ll999l), but similar to the val - 
ues found for most star-forming galaxies (Malhotra et alj|200ll) . 
The peak in the [Ol] 63 ;um/[Cll] 158 fim map coincides with the 
N peak in the PACS 70//m map, suggesting that at this location, 
a star formation rate density higher than in the rest of the ring 
is enhancing both dust emission and [Ol] 63 jim line emission. 
The ratio is also enhanced on a region where the ring and the 
dusty spiral arms intersect, which indicate a possible shock con- 
tributio n to the [Ol] 63 um flux. Mo dels of emission lines from 
shocks (i HoUenbach & McKeJll989h predict values of the [Ol] 
63 /im/[Cll] 158 //m ratio of at least 10. Our measured value of 
~ 0.4 are at odds with a pure-shock interpretation of the line 
ratios. On the other hand, line emission from shocks cannot be 
ruled out because the region in the PACS beam could be a com- 
posite of both shocked gas and additional [Cll] flux contributed 
from star formation. Models of line emissi on from photodissoci- 
ation regions (PDRs) around O stars (e.g.. lKaufman et anil999l) 
show that [Ol] 63 ;um/[Cll] 158 fim increases with the radiation 
field intensity Go and neutral gas density n. However, at low den- 
sities, [Ol] 63 /im/[Cll] 158 /im is mostly sensitive to Go. Using 
iKaufman et al. U999) models for n < 10** cm'"*, we estimat e 
Go to he between 60 - 300, and using iMeiierink et all (l2007h . 
Go is estimated betwe en 200 - 600. The higher Go found using 
IMeiierink et alJ (l2007h is likely due to their lower adopted O/C 
abundance ratio. 

The [Nil] I22ij.m/[NII] 205 A^m ratio varies between 4.0-6.0 
throughout the central 40". The values that correspond to the 
[Nil] 122 ^m peaks are 4.8 and 4.5, whereas at the [Nil] 205 
fim peak, the [Nil] 122 //m/[Nll] 205 //m ratio is 4.0. From the 
comparison with the 70yum contours, we see that the highest val- 
ues of the [Nil] 122 yum/[Nll] 205 fj.m ratio coincide with the 
ring, while the lowest are found in the inner region of the ring 
and also where the dust lanes meet the ring. The variations in 
the [Nil] 122 yum/[Nll] 205 yum ratio are therefore due to a vari- 
ation in the ionized gas density between the ring and the inner 
region. Using a five level model of N^ we find that the variation 
of [Nil] 122 //m/[Nll] 205 fim ratios between the inner region 
and the ring corresponds to a variation of the electron density 
from 150 to 400 cm"-'. The results are insensitive to typical gas 
temperatures (T ~ 6000 - 10000 K) in photoionized gas. This 
means that the density increases by at least a factor of 5 in the 
ring compared to the inner region and the region where the dust 
lane and the ring meet. These values are similar to the central 
region of M82, in which a mean ratio value of 4.2^ [^ was mea- 
sured, corresponding to a mean electron density of 180^j2q cm 
"^ across the central 50" (iPetuchowski et alll 19941) . This is also 
consistent with the value of ~ 220 cm"^ fro m the mid-infrare d 
[Sill] 18/im/33yum ratio over the same region dPale et al.ll2007h . 

The peak in the [Nil] 205 fim emission line map in Fig. 2 has 
a line flux of ~ 0.1 1 x 10"^ W m'^ sr"' and a [Cll] 158 jum/[Nll] 
205 jum ~ 45. For 150 < «e < 400 cm"^ the [Cll] 1 58 ;um/[Nll] 
205 fim ratio in ionized gas is expected to be ~ 3 dOberst et alj 
12006) and thus most (~ 90%) of the [Cll] 158 fim emission we 
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Fig. 5. Top: map of the [Ol] 63 jum/[Cll] 158 fim ratio. The [Ol] 
63 fim map was smoothed to match the resolution at 158yum. 
Bottom: map of the [Nil] 122 fim/[mi] 205 fim ratio. The [Nil] 
122 fim map was smoothed to match the [Nil] 205 fim map. Both 
maps are overlaid with contours from the PACS 70jum map, after 
smoothing to the resolution of each map. All the ratio maps were 
built from line maps clipped at a 2cr level above the noise. The 
cross marks the location of the nucleus. 



measure is coming from neutral gas. The gas heating efliciency, 
measured by the ([Ol] + [Cll]) /FIR ratio, seems to stay constant 
in the mapped area. The log of the ratio is ~ -2.2 inside the ring 
and ~ -2.3 on the ring, both consi stent with the values found in 
nearby galaxies bv lMalhotra et al.l (|2001|) . 

In summary, we have used the PACS Spectrometer to map 
the [Ol] 63 fim, [Olll] 88 fim, [Cll] 158 fim, [Nil] 122 fim and 
206/im far-infrared cooling line emission in the central 5 kpc 
of NGC 1097 for the first time. While the [Ol] 63 fim, [Olll] 
88 fim and [Nil] 122 fim line maps appear qualitatively similar, 
the [Nil] 205 fim map shows a different distribution. The [Ol] 63 
//m/[Cll] 158 fim map shows a relative hotspot on the NE portion 
of the ring indicative of a stronger radiation field or a region 
of shocked gas. The [Nil] 122 yum/[Nll] 205 fim map shows a 
clear increase of ionized gas density in the ring, associated with 
massive star formation activity. 
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